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Abstract 

We survey some of the recent developments in the extraction and application of heavy quark Parton Distribution 
Functions (PDFs). We also highlight some of the key HERA measurements which have contributed to these advances. 
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Two Decades of HERA physics 

The HERA electron-proton collider ring began its 
physics program in 1992 and completed accelerator op- 
erations in 2007. The data collected by the HERA fa- 
cility allowed for physics studies over a tremendously 
expanded kinematic region compared to the previous 
fixed-target experiments. This point is illustrated in Fig- 
ure [T] where we display the e + p Neutral Current (NC) 
cross section vs. Q 2 for the HERA data (runs I and II) 
together with the fixed-target data. We observe that the 
HERA data allows us to extend our reach in Q 2 by more 
than two decades for large to intermediate x values, and 
also extends the small x region down to ~ 10 . 

Additionally, the large statistics and reduced system- 
atics of the experimental data demand that the the- 
oretical predictions keep pace. Over the lifetime of 
HERA we have seen many of the theoretical calcula- 
tions advanced from Leading-Order (LO), to Next-to- 
Leading-Order (NLO), and some even to Next-to-Next- 
to-Leading-Order (NNLO). 

As the required theoretical precision has increased, it 
has been necessary to revisit the many inputs and as- 
sumptions which are used in the calculations. We will 
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Figure 1: e + p NC cross section for the combined HERA data as 
compared with the HERAPDF1.5 fit as a function of Q- for different 
values of X. (Figure from HI prelim- 10- 142 & ZEUS-prel-10-01S.) 
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Figure 2: The computed nuclear correction ratio, F^ e /F^ , as a function of x for Q 2 = 5 GeV 2 . Figure-a) shows the fit (fit B from Ref. 1 16 |) using 
charged-lepton-nucleus (f ± A) and DY data whereas Figure-b) shows the fit using neutrino-nucleus (vA) data (fit A2 from Ref. |2|). Both fits are 
compared with the SLAC/NMC parameterization [T5), as well as fits from Kulagin-Petti (KP) (Ref. and Hirai et al. (HKN07), (Ref. 0). 

The data points displayed in Figure-a) come from a selection of SLAC and BCDMS data. lo1l7ll^l^ [T0l[TT1[T2l . 



examine the role that the heavy quarks-and their associ- 
ated masses-play in these calculations, both for the evo- 
lution of the parton distribution functions (PDFs) and 
also the hard-scattering cross sections. 

Determining the Heavy Quark PDFs 

HERA's reach to larger Q 2 and smaller x values takes 
us to a new kinematic region where the heavy quarks 
(s,c,b) play a more important role. For example, the 
strange and charm quark contributions to Fi at small 
x values can be 30% or more of the total inclusive re- 
sult. To make high-precision predictions for the struc- 
ture functions we must therefore be capable of reducing 
the uncertainty of these heavy quark contributions; this 
requires, in part, precise knowledge of the PDFs which 
enter the calculation. 

The determination of the PDFs requires a variety 
of data sets which constrain different linear combina- 
tions of the PDF flavors. For example, Neutral Cur- 
rent (NC) charged-lepton Deeply Inelastic Scattering 
(DIS) (at low Q 2 ) probes a charge weighted combina- 
tion ~ Au + d + s + Ac. In contrast, charged current 
(CC) neutrino DIS can probe different flavor combina- 
tions via W ± -boson exchange; additionally the neutrino 
measurements can probe the parity-violating xFj, struc- 
ture function. 

Nuclear Correction Factors 

The most precise determination of the strange quark 
PDF component comes from neutrino-nucleon (v-N) 
di-muon DIS (vN — > [T^X) process. This (domi- 
nantly) takes place via the Cabibbo favored partonic 



process vs — » pTc followed by a semi-leptonic charm 
decay. As the neutrino cross section is small, this mea- 
surement is typically made using heavy nuclear targets 
(Fe, Pb), so nuclear corrections must be applied to relate 
the results to a proton or isoscalar nuclei. 

However, recent analyses indicate that the nuclear 
corrections for the v-N and f=-N DIS processes are dif- 
ferent [ 1 1; hence, this introduces an uncertainty into the 
strange quark PDF extraction which was not realized 
previously. Figure[2]displays the nuclear correction fac- 
tors obtained for the a) £*-N and b) v-N processes. Here, 
we plot the ratio of F^ e to an isoscalar F^ as a function 
of x for the Q 2 value indicated. 

The contrast between the charged-lepton (^*A) case 
and the neutrino (vA) case in Figure [2] is striking; 
while the charged-lepton results generally align with 
the SLAC/NMC 0, KP EH) and HKN determi- 
nations, the neutrino results clearly yield different be- 
havior in the intermediate x-region. We emphasize that 
both the charged-lepton and neutrino results are not a 
model — they come directly from global fits to the data. 
To emphasize this point, we have superimposed illus- 
trative data points in the figures; these are simply a) the 
SLAC and BCDMS data llSl 171 l8l l9l IT01 fill Hill orbl 
the vA DIS data [ 13 1 scaled by the appropriate structure 
function, calculated with the proton PDF of Ref. Q. 

The mis-match between the results in charged-lepton 
and neutrino DIS is particularly interesting given that 
there has been a long-standing "tension" between the 
light-target charged-lepton data and the heavy-target 
neutrino data in the historical fits lfl4l fT3Tl . This 
study demonstrates that the tension is not only be- 
tween charged-lepton light-target data and neutrino 
heavy-target data, but we now observe this phenomenon 



K. Kovafik et al. Proceedings of the Ringberg Workshop: New Trends in HERA Physics (2011 ) 



3 




Figure 3: The a) gluon xg(x, Qq) and b) strange quark xs(x, Qq) nuclear PDFs as a function of x for a selection of nuclear A values 
{1,2,4,9, 12,27,56, 108,2071 (from top to bottom at x = 0.01). We choose Q = 1.3 GeV. 



in comparisons between neutrino and charged-lepton 
heavy-target data. 

The nCTEQ PDFs 

The above example underscores the importance of a 
comprehensive treatment of the nuclear corrections to 
achieve the precision demanded by the current preci- 
sion data. To move toward this goal, the nCTEQ project 
was developed to extend the global analysis framework 
of the traditional CTEQ proton PDFs to incorporate a 
broader set of nuclear data thereby extracting the PDFs 
of a nuclear target. In essence, a nuclear PDF not only 
depends on the momentum fraction x and energy scale 
Q, but also on the nuclear "A" value: f(x, Q,A). The 
structure of the nCTEQ analysis is closely modeled on 
that of the proton global analysis; in fact, the nuclear 
parameterizations are designed efficiently to make use 
of the proton limit (A=l) as a "boundary condition" to 
help constrain the fit. 

In Figure|3]we display the gluon and strange nuclear 
PDFs as a function of x for a selection of nuclear A 
values. We observe that for x ^ 0.01 the nuclear modi- 
fications for the strange quark can be ~25%, and for the 
gluon can be even larger. The details of the nuclear PDF 
analysis is discussed in Refs. flT] |T6 2^] 

The nCTEQ web-page contains 19 families of nPDF 
grid files which may be used to explore the variation 
due to the different data sets and kinematic cuts. In par- 
ticular, there is a collections of nPDFs which interpo- 
late between that of Figure[2]-a) which uses the charged- 
lepton-nucleus (PA) data and of Figure[2]-b) which uses 
the neutrino-nucleus (vA) data. 
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2 The nuclear PDFs are available on the web from the nCTEQ page 
at http://projects.hepforge.org/ncteq/ which is hosted 
by the HepForge project. 



Figure 4: k(x) vs. x for 2 = 1-5 GeV for a selection of PDFs, where 
we define k(x) = 2s(x) / (u(x) + d(x)) 



The Strange Quark PDF 

We now compare a selection of s(x) distributions to 
gain a better understanding of the uncertainties arising 
from the nuclear correction factors used to analyze the 
vN DIS . One measure of the strange quark content of the 
proton is to compare s(x) with the average up-quark and 
down-quark sea PDFs: (ii(x) + d(x))/2. Thus, we define 
the ratio k(x) = 2s(x)/(S(x) + d(x)). If we had exact 
SU(3) flavor symmetry, we would expect k = 1; the 
extent to which k is below one measures the suppression 
of the strange quark as compared to the up and down 
sea. In Figure |4] we display k(x) for some recent CTEQ 
PDFs and note that k(x) has a large variation, especially 
at small x values. This reflects, in part, the fact that the 
strange quark is poorly constrained for x < 0.1. For the 
CTEQ6, 6HQ, 6.1, and 6.5 PDF sets, the strange quark 
was arbitrarily set to ~ 1/2 the average of the up and 
down sea-quarks. For the CTEQ6.6 PDF set, the strange 
quark was allowed additional freedom; this is reflected 
in Figure [5] which compares the relative uncertainty of 
the strange quark in the 6. 1 and 6.6 PDF sets. 
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Figure 5: Relative uncertainty of the strange quark PDF as a function 
of x for Q = 2 GeV. The inner band is for the CTEQ6.I PDF set, and 
the outer band is for the CTEQ6.6 PDF set. The band is computed as 
Sj(x)/so(x) where sq(.x) is the central PDF for each set; for CTEQ6.1, 
i = [1,40], and for CTEQ6.6, i = [1,44]. 



W* Production ( y/s = 14 GeV) 
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Figure 6: The differential cross section (dcr/dy) for W + production 
at the LHC (-y's=14TeV) as a function of rapidity y. The partonic 
contributions are also displayed. At y = 0, the contributions (from top 
to bottom) are [total, ud, cs, us, cd } 

W/Z at the LHC 

The above reexamination of the nuclear corrections 
introduces additional uncertainties into the data sets, 
which manifests itself in increased uncertainties on the 
strange quark PDF. To see how these uncertainties might 
affect other processes, we consider, as an example, W/Z 
production at the LHC. As we go to higher energies, the 
heavy quarks will play an increasingly important role 
because we can probe the PDFs at smaller x and larger 
Q; this means that the heavy quark PDF uncertainties 
can have an increased influence on LHC observables 
compared to Tevatron observables. 

In Figure [6] we display the LO differential cross sec- 
tion for W production at the LHC as a function of rapid- 
ity y, as well as the individual partonic contributions. 
We note that in the central rapidity region the contribu- 
tion from the heavy quarks can be 30% or more of the 
total cross section; this is in sharp contrast to the sit- 



uation at the Tevatron where the heavy quark contribu- 
tions are minimal. Thus, a large uncertainty in the heavy 
quark PDFs can influence such "benchmark" processes 
as W/Z production at the LHC. Of course, given the high 
statistics from the LHC (the 201 1 proton-proton run ex- 
ceeded 5 fb~ l ), it may be possible to turn the question 
around and ask to what extent the LHC data may con- 
strain the heavy quark PDFs. 

Zero Mass (ZM) and General Mass (GM) Schemes 

We now turn to charm production and the measure- 
ment of the charm PDF. HERA extracted precise mea- 
surements of F c 2 and F*, and recently these analyses 
have been updatecj^jto include the low Q 2 data to cover 
the kinematic range of Q 2 = [2, 1000] GeV 2 and x down 
to 10~ 5 . E1ED 

A global fit of HERA I data flS ED for F c 2 was per- 
formed using both the General Mass Variable Flavor 
Scheme (GM-VFS), and also the Zero Mass Variable 
Flavor Scheme (ZM-VFS). (2T] While the GM-VFN re- 
sult yielded an improved^- 2 , the ZM-VFN results-when 
implemented consistently-yielded an acceptable fit to 
the data. Given the expanded kinematic coverage of the 
recent HERA data, it would be of interest to repeat this 
comparison. Presumably the new data sets would allow 
for increased differentiation between the ZM-VFN and 
GM-VFN scheme results. 

Choice of Theoretical Schemes 

Having illustrated the impact of different theoretical 
schemes on the data analysis, we take a moment to com- 
pare and contrast some of the different schemes that are 
currently being used for various PDF analysis efforts. 
While many of the global analyses use a Variable Flavor 
Number (VFN) scheme to include the heavy quark as a 
parton, the detailed implementation of this scheme can 
lead to notable differences. At the 2009 Les Houches 
workshop, a comparison was performed among a num- 
ber of the different programs to quantify these differ- 
ences. All programs used the same PDFs and as values 
so that the differences would only reflect the particular 
scheme. The complete details can be found in Ref. [22 1, 
and Figures [7] and [8] display sample comparisons. 

Figure[7]compares the S-ACOT scheme which is used 
for the CTEQ series of global analyses, 11231 l24l and 
the FONLL which is used by the Neural Network PDF 
(NNPDF) collaboration. [25 1 In the figure, these two 



3 Cf., ZEUS-prel-09-015 
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Figure 7: Comparison of F% for the Fixed-Order-Next-to-Leading- 
Log (FONLL) with the Simplified-ACOT (S-ACOT) scheme. There 
are four curves displayed, two for the ordinary ^-rescaling, and two 
for an alternate ^-rescaling (labeled "v2"). The FONLL and S-ACOT 
results are identical throughout the x range. (Figure from Ref. i22t .j 
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Figure 8: Comparison of F% for the Fixed-Order-Next-to-Leading- 
Log (FONLL) and the MSTW08 NLO results. The FONLL results 
are shown for both the "A" and "B" variations. The FONLL results 
differ from the MSTW08 results for low Q and x; for larger values of 
Q and x they are more comparable. (Figure from Ref. K22\ . ) 



implementations (S-ACOT and FONLL-A) are numeri- 
cally equivalent. 

Figure [8] compares two variations of the FONLL 
scheme ("A" and "B") with the MSTW08 results which 
is used in the MSTW series of PDF global analy- 
ses. Il26l l27ll Here these differences reflect the differ- 
ent organization and truncation of the perturbation ex- 
pansion; it does not indicate that one choice is right or 
wrong. We expect such differences to be proportional 
to ~ a N s x 0(m 2 /Q 2 ). Thus, as we increase the order of 
perturbation theory or the energy scale the differences 
should decrease; we have explicitly verified the differ- 
ence is reduced as Q 2 increases, as it should. When we 
are able to carry these calculations out to higher orders, 
the scheme differences should be further reduced; this 
work is in progress. 

Charm Mass Dependence and F c 

The experimental extraction of the "inclusive" F c 2 re- 
quires a differential NLO calculation of DIS charm pro- 
duction to be extrapolated over the unobserved kine- 
matic regions. These analyses generally make use of the 
HVQDIS program 11281 which computes F c 2 in a Fixed- 
Flavor-Number (FFN) scheme at NLO. In this calcula- 
tion, the charm is produced only via a gluon splitting, 
g — > cc, and there is no charm PDF. Thus, the charm 
mass (m c ) enters only the partonic cross section &(m c ) 
and the final state phase space; there is no PDF charm 
threshold. 

Although we would also like to perform the ex- 
traction of F c 2 using a Variable-Flavor-Number (VFN) 
scheme, the challenge is that no NLO differential pro- 
gram exists for this process. In lieu of a VFN extraction, 
another avenue is to study the influence of different the- 
oretical schemes and m c parameters in the analysis of 
the F c 2 data. We describe such a study below. 

In many analyses, the value of the charm mass is 
taken as an external fixed parameter. A recent inves- 
tigation has taken a closer look at the role of the charm 
mass parameter m c and examined the combined effects 
of m c and the theoretical scheme used; preliminary re- 
sults of this study are displayed in Figure [9] For each 
of the schemes listed in the legend, fits were generated 
for fixed m c values in the range [1.2, 1.8] GeV. Thus, 
the minimum of the x 2 curve represents the "optimal" 
choice of the charm mass parameter m c for that specific 
scheme. 

We observe that the various schemes prefer m c val- 
ues ranging from 1.2 to 1.7 GeV, The largest m c value 
(1.68 GeV) comes from the Zero Mass VFN Scheme 
(ZM VFNS) which only uses m c for the PDF charm 
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Figure 9: Comparison of x for HERA I + FS E fits using different 
heavy flavor schemes as a function of the charm quark mass parameter 
m modd ^ Figure from HI prelim- 10- 143 &■ ZEUS-prel-10-019) 



threshold; it is absent in the phase space for the zero- 
mass case. In contrast, for the S-ACOT-;^ scheme the 
notation^] indicates there are effectively two factors 
of m c in the final phase space, and this yields the small- 
est value of m c (1.26 GeV). 
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Figure 10: Measurement of y + c and y + b vs. Pj as measured by 
D-Zero. (Figure from Ref. [321.) 
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Figure 1 1 : The integrated momentum fraction (in percent) of the "ex- 
trinsic" charm and bottom quarks generated by gluon splitting as a 
function of the scale Q. Reference lines are indicated at 1% and 2%. 



Extrinsic & Intrinsic Charm PDFs 



The ACOT scheme and the Roberts-Thorne (RT) 
scheme yield m c values in the intermediate region. 
The ACOT scheme uses the full kinematic mass re- 
lations in the partonic relations, and the scaling vari- 
able is intermediate (m c =1.58 GeV) between the ZM- 
VFN scheme and the S-ACOT-^ scheme. The S-ACOT 
scheme (not shown) is virtually identical to the full 
ACOT scheme, and also yields a m c value in the in- 
termediate region. Therefore, comparing the S-ACOT- 
X and S-ACOT schemes, we find that the ^-rescaling 
variable is dominantly responsible for the shift of the 
optimal m c value from ~ 1.26 to ~ 1.58. This obser- 
vation suggests that it is the rescaling of the x variable 
which enters the PDFs that generates the dominant ef- 
fect. While this study is continuing, it does indicate 
the sensitivity of the charm mass and scheme choice in 
these precision analyses. 



The charm quark and bottom quark PDFs can be 
probed directly at the Tevatron by studying photon- 
heavy quark final states which occur via the sub-process 
gQ ~ * jQ at LO. This process has been measured at the 
Tevatron for both charm and bottom final states, and we 
display the results in Figure 



10 



as a function of P y T for 
two rapidity configurations. This measurement is par- 
ticularly interesting as the dominant process involves 
a heavy quark PDF; this is in contrast to DIS charm 
or bottom production, for example, where over much 
of the kinematic range the process is dominated by the 
gluon-initiated process (e.g., yg — > QQ) rather than the 
heavy quark initiated process (yQ — > Q). 



4 Specifically, the ^-prescription rescales the partonic momentum 
fraction via x — > x(l + (2m c /Q) 2 ) in contrast to the traditional Bar- 
nett 1291 "slow-rescaling" which is x — » x(l + (m c j Q) 2 ). 
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Figure 14: Measurement of Ft. using the combined HERA data set 
from HI and ZEUS. The data are compared with a selection of the- 
oretical predictions. (Figure HI prelim- 10-044 & ZEUS-prel- 10-008). 



Figure 12: dcr/dp Tr (pb/GeV) at the Tevatron for the CTEQ6.6 PDFs, 
and two intrinsic charm (IC) models. (Figure from Ref. HOf . ) 
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Figure 13: Ratio of dcr/dp Ty (pb/GeV) at the LHC ( V* = 1 TeV) for 
the BPHS IC model to the CTEQ6.6 PDF for a selection of rapidity 
bins. ( Figure from Ref. tfOV . ) 



Examining Figure 10 we observe that the bot- 



tom quark production measurements compare favor- 
ably with the theoretical predictions throughout the P y T 
range, but the charm results rise above the theory pre- 
dictions for large P 7 T . Although there may be a num- 
ber of explanations for the excess charm cross section 
at large P y T , one possibility is the presence of intrinsic 
charm (IC) in the proton. In the usual DGLAP evolu- 
tion of the proton PDFs, we begin the evolution at a low 
energy scale Qq < m c and evolve up to higher scales. 
The charm and bottom PDFs are defined to be zero for 
Q < m c j,, and above the mass scale the heavy quark 
PDFs are generated by gluon splitting, g — > QQ; we 
refer to this as the "extrinsic" contribution to the heavy 



mentum fractions for the charm and bottom quarks as a 
function of the scale Q; these are zero for Q < m c j,, and 
then begin to grow via the g — > QQ process. 

It has been suggested that there may also be an "in- 
trinsic contribution to the heavy quark PDFs which is 
present even at low scales Q < m Cy b- While it is difficult 
to constrain the detailed functional shape of any intrin- 
sic heavy quark distribution, the total momentum frac- 
tion of any intrinsic contribution must be less than ap- 
proximately 1% if it is to be compatible with the global 
analyses. 



In Figure 12 we illustrate the effect of including an 
additional intrinsic charm component in the proton. The 
BHPS IC model concentrates the momentum fraction at 
large x values, and the Sea-like IC model distributes the 
charm more uniformly^] It is intriguing that the IC mod- 
ification of the proton PDF can increase the theoretical 
prediction in the large P y T region, but this observation 
alone is not sufficient to claim the presence of IC; this 
would require independent verification. In Figure 13 
we display the cross section ratio for y + c at the LHC 
for the BHPS IC model for a selection of rapidity bins. 
Thus the LHC can validate or refute this possibility with 
a high-statistics measurement of y + c, especially if they 
can observe this in the forward rapidity region. 

The Longitudinal Structure Function F L 

Most of the previous discussion has addressed the de- 
termination of the quark PDFs. Constraining the gluon 
PDF is a challenge, and the longitudinal structure func- 
tion Fi is particularly interesting as it involves both the 



quark PDFs. In Figure 1 1 we display the integrated mo- 



s For details, c.f., Refs. 13011311 
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a) 



b) 



Figure 15: Fractional flavor decomposition of F'JFl, vs. Q in GeV for a) x — 10 and b) x = 10 . Reading from the bottom, we plot the 
cumulative contributions for [u, d, s, c, b). 



heavy quark and the gluon distributions. Using the com- 
bined data from HI and ZEUS, HERA has extracted Ft, 
in an extended kinematic regime, and Figure[T4]displays 
the result of the combined data as compared with vari- 
ous theoretical predictions. 

The measurement of F L is special for a number of 
reasons, and we write this schematically as: 



Ft. - 



Q- 



q(x) + a s [C g ® g(x) + C q ® q(x)} . ( 1 ) 



Note that the LO term is zero in the limit of massless 
quarks as the (m 2 /Q 2 ) factor in Eq. ([T]i suppresses the 
helicity violating contributions; this is a consequence of 
the Callan-Gross relation. Therefore, for light quarks 
the dominant contributions come from the NLO gluon 
term; hence, Ft, can provide useful information about 
the gluon PDF. 

For the heavy quarks the picture is less obvious. 
While the NLO heavy quark contributions will clearly 
be small compared to the dominant gluon terms, the 
heavy quarks can contribute at LO if they can overcome 
the (m 2 /Q 2 ) suppression. This is why the prediction of 
Fi into the low Q 2 region as measured in Figure 14 is 
such a theoretical challenge. This raises a number of 
questions: What is the flavor composition of FjP. Where 
are the heavy quark contributions important? 



In Figure 15 we display the fractional contributions 
to the structure functions F'JF L vs. Q. We observe that 
for large x and low Q the heavy flavor contributions are 
minimal. For example, in Figure 15 a) at Q ~ 5 GeV we 
see the M-quark structure function F U L comprises ~ 80% 
of the total, F d L is about 10%, and the s, c and b quarks 
divide the remaining fraction. 

At smaller x values the picture changes and the heavy 
quarks are more prominent. In Figure 15 b) for Q ~ 



2 GeV we see the M-quark structure function F" L com- 
prises ~ 55%, F d L and F^ are both about 20%, and the 
c and b quarks make up the small remaining fraction. 



However, F° L increases quickly as Q increases and is 
comparable to F'[ (~ 40%) for Q ~ 20 GeV. Addition- 
ally, for large Q ~ 100 GeV we see the contributions 
of the M-quark and c-quark are comparable, the li-quark 
and s -quark are comparable, and the relative sizes of the 
u, c to d, s terms are proportional to their couplings: 4/9 
to 1/9. Thus, for low x and intermediate to large Q val- 
ues we see that the quark masses (aside from the top) no 
longer play a prominent role and we approach the limit 
of "flavor democracy." 

Concluding Remarks 

We reviewed a number of recent developments re- 
garding the extraction and application of heavy quark 
Parton Distribution Functions (PDFs). The high pre- 
cision HERA measurements were essential in develop- 
ing and refining the theoretical treatment of the heavy 
quarks. Even though the accelerator facility stopped op- 
eration four years ago, the analysis of the data contin- 
ues. The results of these analyses will provide the foun- 
dation upon which future PDF analyses will be built, 
and the advances of the experimental analysis and theo- 
retical tools developed at HERA will continue to influ- 
ence future hadronic studies including those now begin- 
ning at the LHC. 
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